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Abstract. We derive a Khinchine-Pollaczek formula for random 
walks whose steps have a geometric left tail. The construction 
rests on the memory-less property of the geometric distribution. 
An example from a tandem queue modeling dynamic instability in 
microtubules is given. 



1. Introduction 

Many results in queueing theory are based on the connection of the 
waiting times process and an associated random walk. If the queue 
is an M/M/1 queue, then a classical result for the associated random 
walk is the Khint chine- PoUaczek formula, [4, XI 1. 5]. The derivation 
exploits the "memoryless" property of the exponential distribution. 
In this paper we obtain a discrete version of this result, using the 
memoryless property of the geometric distribution. The discreteness 
of the distribution introduces a number of additional subtleties, but in 
general our proof is organized as the proof of the continuous case in [2] . 

The motivation to consider the discrete case came from our attempt 
to model what is known as "dynamic instability" of microtubules (see 
PQ Ch. 16]) by a tandem queue. While the primary queue of the 
tandem is M/M/1, the secondary queue is not. However, the primary 
queue bequeathes a geometric left tail to the secondary queue and this 
caused us to search for a Khintchine-PoUaczek formula in that setting. 



2. Preliminaries 

Let Xi, X2, . . . , be a sequence of independent random variables with 
common distribution F not concentrated on a half-axis. The induced 
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random walk is the sequence of random variables 

5*0 = 0, Sn = Xi + ■ ■ ■ + Xn- 

The integer n is a (strict) ladder index for the random walk if 

max Sk < Sn- 

0<k<n 

If n is a ladder index, then Sn is the ladder height associated with n. 
For Borel sets A C (0, oo), define a finite measure L by 

(1) L(A) = V P f max Sk = < Sn E a] . 

n=\ ^ - / 

The probability that there is at least one ladder index is 

(2) p = L(0,oo) = P (sup5„ > 



Let Ti be the first ladder index, and denote ifi the first ladder height, 
i.e. li\ = S'tj. These variables are defective with probability 1—p and 
remain undefined if there is no first ladder index. In fact, for x > 0, 

(3) L[x,oo) = P{Hi>x). 

Following the notation in ^ Chapter XII], we call the smallest n 
such that Si < 0, ... , Sn-i < 0, but S'„ > the first weak ladder index 
and denote it by Ti. The corresponding weak ladder height is denoted 
by Hi, so that Hi = Sj^^. Again, these variables are possibly defective. 
We set 



(4) 



( = p[Hi = 0) = ^P ( max Sk < 0, 5„ = o\ 

n=l \1^*^<"' / 



Since Xi > implies Hi > and it is assumed that Xi is not concen- 
trated on a half-axis, it follows that < C < 1. 

Lemma 1. For a; > 0, we have 

(5) P (Ti < Ti,Hi >x)= CP{Hi > x) . 

Proof. If < /c < n and P (Ti = k, Hi = O) > 0, then 

P (Ti = n,Hi> x\Ti = k,Hi = O) = P {Ti = n - k, Hi > x) 
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by the Markov property. Thus 
P{Ti<T,,H,>x) 



oo n—1 



= J] J]P (Ti = A;,Ti = > x) 

n=2 k=l 

oo oo 

= J2 Yl ^ (^1 = ^'^1 = 0) P (Ti = n,Hi> x\Ti = k,Hi = O) 

k=l n=k+l 

(6) 

oo oo 
k=l n=l 

□ 

Denote L"* the n-fold convolution of L with itself and a unit 
mass at the point 0. Define the measure ijj by 

oo 

(7) ^(A) = 5^L"*(A), Ac[0,oo). 

n=0 

We have, cf. p, Theorem 24.2(iii)], 

Theorem 2. //p < 1, then with probability p"'{l —p) there are exactly 
n ladder indices; with probability 1 there are only finitely many ladder 
indices and sup„ Sn < oo; finally 

(8) P(sup5„ G A) = (l-p)^(A), Ac[0,oo). 

Vn>0 / 

Furthermore, cf. [21 Theorem 24.3], 
Theorem 3. The measure if) satisfies 



(9) 



f ^[0,a;-y] dF{y) = Y P ( mm Sk>^,Sn<x\ . 



3. One-sided geometric tails 



We assume from now on that the random variables Xi are integer 
valued. 

Theorem 4 (Geometric right tail). Suppose that E[Xi] < and that 
the right tail of F is geometric: 

(10) P{X^>x)=ir\ a; = 0,1,..., 
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where < ^ < 1 and < r < 1. Then p < 1 and 

(11) P (sup Sn>x] ^p[l-{l-p){l- r)]" , X = 0, 1, ... . 

\n>0 / 

Moreover, 1/[1 — (1 — p)(l — r)] is the unique root of the equation 

oo 

(12) Yl ^"^(^1 = x) = l 

x=—oo 

in the range 1 < s < 1/r . 
Proof. For x > we have 

P I max Sk < 0, Sn > x ] 

\^0<fc<n J 

(13) ^ p( Sn> xlSn > 0, max Sk<o]p ( max Sk < 0, > ) , 

\ 0<A;<n J \0<k<n J 

and also 



P { Sn> x\Sn > 0, max Sk<0 

\ 0<k<n 

(14) = P ( X„ > X - > -Sn-u max 5fc < ) . 

Furthermore, since Xn is independent of Sq, . . . , Sn-i, and the geomet- 
ric distribution is memoryless, 

P {Xn> X- Sn-l\Xn > Sn-l, maX S'fe < 
y 0<fc<n 

(15) = ^P{Xn >x)^ r\ 
Whence 

(16) P ( max Sk<Q,Sn> x] ^r'^P ( max -5^ < < 5„ | , 

\0<ifc<n J \0<k<n J 

and so 

L[{x, oo)) = y PI max S'fe < 0, 5'„ > x ) 

n=l ^ ~ ^ 

(17) = p f max 5fc < < = r^p. 

n=l ^ ~ / 
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This shows that the measure L is p times the geometric distribution 
with parameter 1 — r. Hence L"* is times the negative binomial 
distribution with parameters n and 1 — r: For x > 0, 



{lZ\) (1 - r)"r^-", if 1 < n < 
0, otherwise. 



L"*{a;} 

In particular, for x = 1, 2, . . . , 

(1 - r)V^-" = p{l - r)[l - (1 - p){l - r)]^-\ 
SO that 

X 

^([0, x]) = 1 + Y,p{l - r)[l - (1 - - r)]^-i 



ra=l ^ 



:i8) 



fc=i 

1 + (1 — r)x, if p = 1; 

T^-T^[l-(l-|')(l-Or, ifp<l- 



Since E[Xi] < 0, we have Sn — — oo P-a.s and p < 1. The identity 
(fTTj) now follows from ([8]). 

Next, since Sn — oo, the right side of ([9]) is equal to one for x = 0, 
so that by ([18]) 



(19) (1 - P [1 - (1 - - r)r) F{y} = l-p. 

y=-oo 

Denote f{s) the probability generating function of Xi given on the left 
in f|T2|) . Because of ffTOl) . it is defined for < s < 1/r and 

^-^ 1 — rs 

j/=-oo 

Since 1 < 1/[1 — (1 — p)(l — r)] < 1/r, it follows from this, (fT9|) . and 
F(-l) = 1 - that /(1/[1 - (1 - p)(l - r)]) = 1. 
Finally, /(I) = 1 and, since P(Xi < 0) > 0, 

/"(s) = E x(x - l)s''-^F{x} > 

Z\{0,1} 

for < s < 1/r. Hence equation (fT2l) cannot have more than one root 
in (1, 1/r). 

□ 
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Theorem 5 (Geometric left tail). Suppose that E[Xi] < and that the 
left tail of F is geometric: 



(20) P{Xi <x)= ^r"^, X = 0, -1, ... , 
where < ^ < 1 and < r < 1 . Then 

(21) (l-C)p = r + (l-r)E[Xi], 
and for x = 1,2, . . . , 

X 

(22) (1 - C)i:(0, x] = F{x) - F(0) + (1 - r) ^(1 - F{m)). 

m=l 

Proof. Assume first not ( 120|) but ( !T0|) . and assume also that E[Xi] > 0. 
In that case, Sn oo, P-a.s., so that p = 1 by ([2]). Thus, from ( !T8|l . 
we have '?/'[0, a^] = 1 + (1 — r)x for a; > 0. Whence, for x > (cf. ([9])), 

EP I min Sk > 0, Sn < X ] 
\l<k<n J 
n=l ^ - ' 

= ^ [1 + (1 - r){x - m)\ P{Xi = m) 

(23) =[! + (!- r)x]P(Xi < x) - (1 - r) ^ mP(Xi = m). 



m<a: 



If (1201) holds with E[Xi] < 0, then the above applies to the sequence 
{—Xn}- In particular, if X^ is replaced by — x by — x, then 
becomes 

EP I max S'fc < 0, S'„ > X I 

n=l ^ - / 

(24) = [1 - (1 - r)x]P(Xi > x) + (1 - r) ^ mP(Xi = m). 



m>a: 



Recall the definitions of the weak/strict first ladder index. We have 
for X > 



P max Sk = < Sn, •S'n > x — P max 5*^ < 0, Sn > x 

\0<k<n J \l<k<n 



(25) =p[Ti<Ti = n, Sn>x). 
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Thus, from ([5]) and ([3]), 

L[x, oo) = ^ P ( max 

n=l V^-'^^" / 

(26) 

= > PI max Sk < ^-1 Sn> x] + C-^f^;? oo)- 

\ l<k<n I 
n=l ^ - ^ 

Since 

^ mP(Xi = m) = a;P(Xi > x) + ^(1 - P(m)), 

ra>^x m>a; 

it follows from (EH) and ([26]) that 

(27) (1 - QL{x, oo) = 1 - P(x - 1) + (1 - r) J](l - P(m)). 

ra>x 

From (EDI), 



(28) ^ mP(Xi = m) = -^r(l - r) J] mr'""^ = -i- 



1 — r 

m<0 m>l 

Using P(0) = this gives 
(29) 

^(1 - P(m)) = E[Xi] - ^ mP(Xi = m) = E[Xi] + ^— P(0). 

m>0 m<0 

Hence 

(1 - C)i^[l, oo) = 1 - P(0) + (1 - r) 5^(1 - P(m)) 

(30) =r + (l-r)E[Xi], 

from which fl2T]) follows. Applying fl30|) and fl27l) to 

L(0, x] = oo) — L\x + 1, oo) 

gives ([22D. □ 

Corollary 6 (Khinchine-Pollaczek formula). Under the assumptions 
of Theorem\^ if 

M(s) = E [s™P">o5"] , = E[s^i,Xi > 0], 

then 

1 - C - r - (1 - r)E[Xi] 



(31) M{s) 



i-C-[(i-le^)-^+(^) + ^i5^(i-i^(o))]- 
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Proof. It follows from ((22]) that 

(1 - OL{x} = P{Xi = x)+ rP{Xi > x). 
Thus, if C{s) = E[s-^^, sup„>o Sn > 0], then 

oo oo oo 

{l-0C{s) = Y,s^P{X, = x) + {l-r)J2s^ J2 P{X, = m) 

x=l x=l m=x+l 

oo 

= J'+(s) + (1 - r) ^ P(Xi = m) 

m=l 

(32) = (l - ^+(.) + .i^(l - F(0)). 
From ([8]) and ([7]) we see that 

oo 

(33) A^(.) = (l-p)5^£"(s). 

71=0 

Using ([32]) and (EID, the result follows. 

□ 

4. An example from queuing theory 

We now apply the above results to a particular distribution F with 
a geometric left tail. The choice of F is motivated by an application of 
queueing theory to the dynamics of cellular filaments. 

4.1. Nucleotide hydrolysis and cap dynamics: a M/M/l-queue. 

ATP (GTP)-bound monomers arrive at the + end of an actin or tubu- 
lin filament and attach themselves. The filament acts as an ATPase (or 
GTPase) and hydrolyses the bound ATP to ADP (or GTP to GDP). 
We assume that hydrolysis is "vectorial," meaning that monomers are 
hydrolyzed one after the other and in the order they arrived. Conse- 
quently, a filament consists of a hydrolyzed part which includes the mi- 
nus end of the filament, and an unhydrolyzed part containing the -I— end 
of the filament. The unhydrolyzed part is called the ATP(GTP)-cap 
(or simply cap). 

The nucleotide-bound monomer at the minus-end of the cap is next 
to a hydrolyzed monomer. We call this transition the hydrolysis edge. 
Only the nucleotide-bound monomer at the hydrolysis edge undergoes 
hydrolysis at that time. We model the dynamics of the cap as a sin- 
gle server queue where the nucleotide-bound monomers are the "cus- 
tomers," the ATP(GTP)-cap is the queue, and the service provided by 
the "server" at the minus end of the cap is hydrolysis. It is assumed 
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that monomers arrive individually at times i„ and that the interar- 
rival times cr„ are independent and identically distributed exponential 
random variables with parameter l/a, so that 

P{(Tn >t} = e"*/", for t > 0, and E[o-„] = a, n = 1, 2, . . . . 

The service time is the duration of hydrolysis (from the moment 
the hydrolysis edge reaches the monomer and hydrolysis commences, 
until hydrolysis of the monomer is completed) for the nth monomer 
arriving after t = 0. The random variables Tn,n = 1,2, ... , are as- 
sumed to be independent, identically distributed exponential random 
variables with parameter and mean E[t„] = /3. Furthermore, we 
assume that the families {<t„, n > 1} and {t„, n > 1} are independent. 

The actual waiting time Wn, n = 1, 2, . . . is the time from the mo- 
ment of arrival of the nth monomer until its hydrolysis commences. 
A busy period is the lifetime of an ATP/GTP-cap, and the number 
of customers present in the system at time t, Xt, is the length of the 
ATP/GTP-cap at time t. Hence {xt = 0} stands for: "at time t the 
filament is uncapped." For an introduction to queueing theory, see [3]. 

4.2. Catastrophe and rescue: a coupled queue. Once filamen- 
tous monomers have undergone hydrolysis, they enter a second queue, 
the hydrolyzed segment of the filament. The nth customer (hydrolyzed 
monomer) joining this queue after t = arrives at time r„, the depar- 
ture time from the first queue (the cap). For this second queue, we 
assume that the "server" is only active when no cap is present. The 
"service" consists of dissociation of a monomer from the filament. Ser- 
vice is provided in the inverse order of arrival — "last come, first served." 
This service policy means that at the moment the ATP/GTP-cap dis- 
appears, hydrolyzed monomers start dissociating from the +-end. Dis- 
sociation continues until an ATP(GTP)-bound monomer arrives and 
"caps" the filament. The service times (dissociation times) Un for this 
second queue are assumed to be independent and identically distributed 
exponential random variables with parameter I/7. The parameter (or 
rate) I/7 is assumed to be much larger than the rate of arrival of 
ATP(GTP)-bound monomers, 1/a, so that when the filament is un- 
capped, it quickly loses monomers — undergoes a "catastrophe" — until 
it is capped again by an ATP (G TP) -bound monomer — the "rescue." 
For an introduction to the dynamics of cellular filaments, see [ll Ch. 
16]. 

Thus a rescue occurs at the beginning of a busy period of the first 
server, while a catastrophe commences when the busy period ends. 
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Denote n the number of monomers added to the second queue during 
a rescue. 

4.3. Results for the ATP(GTP)-cap. Since the cap is modeled as 
an M/M/l-queue, results are readily available in the literature, [SJ II. 2]. 
We assume from now on that a = /3/q; < 1, i.e. that the rate at which 
monomers are hydrolyzed exceeds the rate at which monomers arrive 
at the +-end of the filament. 

• Cap length. The cap length Xt is an ergodic Markov process 
with stationary distribution 

(34) Uj = lira P{xt = j\xQ = k} = lim P{xt = j} = a^{l — a), 
for j = 0, 1, . . . . Thus 



(35) lim E[xt] = , lim Var[xt] = + . 

t^oo 1 - a t^oo (1 - a)2 1 - a 

• Cap lifetime — duration of rescue (busy period). The cap lifetime 
p has mean 

(36) /3/(l-a). 

• Monomers hydrolyzed while filament capped. The probability 
that while the filament is capped exactly k monomers are hy- 
drolyzed is given by 

1 /2A; - 1\ a'^-i 

(37) P{„ = i.} = _( ^ jjrr^' ' = 

These probabilities have the generating function 



,38) t/w^i:p{.^*}.-i±^(i-^i-^ 

which converges for \s\ < 1 + (1 — a)^/4a. It follows readily that 

1 



(39) E[n] = 

• Maximal cap length. The probability that during a busy period 
the maximal cap length does not exceed h monomers is given 
by 
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4.4. The hydrolyzed filament. During a busy period monomers are 
hydro lyzed (serviced in the M/M/l-queue) and become part of the 
hydrolyzed part of the filament ("arrive at" the coupled queue). During 
an idle period of the M/M/l-queue (when the filament is uncapped), 
hydrolyzed monomers dissociate from the +- end of the filament. For 
k = 1,2,..., denote rik the number of monomers hydrolyzed during 
the k-th busy period. The variables nk,k = 1, 2 . . . are independent 
and identically distributed with distribution given by fl37p . 

Since the time until the arrival of the next ATP(GTP)-bound monomer 
is exponential with parameter 1/a and the time until dissociation of a 
hydrolyzed monomer is exponential with parameter I/7, we have 

Pjmonomer dissociates before capping monomer arrives} = . 

a + 7 

Furthermore, since exponential random variables are memoryless, this 
implies that 

P{m monomers dissociate before filament is capped} 

^ m = 0,l,..., 

a + 7/ a + 7 

assuming that the filament contains at least m monomers. 

Denote Zj., k = 1,2, . . . , the length of the filament at the time of the 
k-th rescue. If at t = the length of the filament is 0, then 21 is the 
number of monomers hydrolyzed during the first busy period minus 
the number of monomers dissociated during the first idle period, i.e. 

Zi = max(ni — mi, 0), 

where mi is independent of rii and has distribution given by fj^Tj) . 
More generally, it is easily seen that 

(42) Zk+i = max{zk + rik+i - mk+i,0), k = 0,l,..., 

where the random variables m^, k = 1,2,..., are independent and 
identically distributed with the common distribution ( HTl) . and inde- 
pendent of the variables n^, k = 1,2.... Thus, from [4, VI. 9], we see 
that the sequence Zk,k = 1,2..., can be considered the sequence of 
waiting times for a single server queue, with interarrival times m^, and 
service times n^, k = 1,2, . . . . Since neither variable is exponentially 
distributed, this queue is a G/G/l-queue, for which far fewer explicit 
results of the type listed above for the M/M/l-queue exist. However, 
the interarrival times m^ have a geometric distribution, which is the 
discrete analogue of the exponential distribution. In particular, the 
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interarrival time distribution is also memoryless, 

P{mk > m + n\mk > m} = P{mk > n}, 

for m, n = 0, 1, . . . , and k = 1,2 . . . . 

Since Efn^] = a/{a — 13) (cf. fl39|) ). and E[mjt] = 0/7, it follows that 
the G/G/l-queue has load 

b = ^ 
a — (3 

We have 6 < 1 if and only if /3 + 7 < a, i.e. if the average time for 
hydrolysis and dissociation of a monomer is less than the average time 
between the arrival of ATP(GTP)-bound monomers. 
We now apply our random walk results from above. Let 

Oi 

Xk = Uk- ruk, and r = — ■ — . 

a + 7 

Then 6 < 1 implies 

1 r 

(43) < 0. 

1 — a 1 — r 

Furthermore, 

P(mi >k)= r^. 
Thus, if V denotes the generating function of mi, then 

(44) V{s) = E[s"^i] = \s\ < 1/r. 

Denote F the common distribution of X^. Then F{x) = ^r~^, x = 
0,-1,..., with ^ = U{r), 

E[Xi] = E[ni] - E[mi] = - < 0, 

1 — a I — r 



and 

Moreover, 
and 



r + (1 - r)E[Xi] ^ 



1 — a 

J^+{s) = V{l/s)[U{s)-U{r)], 
1 -C- 

' 1— a 



Mis 



1 - C - (1 - r)si^^ 



A series representation for ( is readily found using Lemma 2 in [U 
XVIII.3] and the fact that the constant term in [U {s)V {1 / s)]"' equals 
P(Xi + --- + X„ = 0). 
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